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Abstract 

We study a series of one-dimensional discrete-time quantum-walk models labeled by half integers 
j = 1/2, 1, 3/2, • • •, introduced by Miyazaki et al, each of which the walker's wave function has 2j+l 
components and hopping range at each time step is 2j. In long-time limit the density functions of 
pseudovelocity-distributions are generally given by superposition of appropriately scaled Konno's 
density function. Since Konno's density function has a finite open support and it diverges at the 
boundaries of support, limit distribution of pseudovelocities in the (2j + l)-component model can 
have 2j + 1 pikes, when 2j + 1 is even. When j becomes very large, however, we found that these 
pikes vanish and a universal and monotone convex structure appears around the origin in limit 
distributions. We discuss a possible route from quantum walks to classical diffusion associated 
with the j — ► oo limit. 
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I. INTRODUCTION 



In contrast with diffusive behavior of classical random walks, in which the standard 
deviation of walker's positions from the starting point is proportional to the square root of 
time-step in the long-time limit, quantum walkers have velocities 

and probability distributions of pseudovelocity, the position divided by time, are discussed 
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In the present paper we consider a series of discrete-time quantum-walk models on the 

n 

one-dimensional lattice (i.e. integers) Z = {• ■ ■ , —2, —1, 0, 1, 2, ■ ■ •}, introduced by [141 ]. 
where models are labeled by half integers j = 1/2, 1, 3/2, • • •. In the model indexed by j, 
the state of quantum walker is described by a (2j + l)-component wave function 



\ 



,(i) 
-in 



^> +1 (x,t) 



x e z, t = 0,1,2,... 



normalized as E xe z \^ j \x,t)\ 2 = 1, where \^( j \x, t)\ 2 



[V®(x,t)]H®(x,t) 



m=—j I Y m 



X 



t)\ 2 . Let be a quantum coin represented by a (2j + 1) x (2j + 1) unitary 



matrix, whose (m, m')-component is denoted by Rl^ m /. In the present paper, when we write 
matrices and vectors whose elements are labeled by m,m', we will assume that the indices 
m and m' run from j to — j in step of —1. At each time step, 2j + 1 components of wave 
function is mixed by the quantum-coin matrix Ry\ and then the quantum walker hops to 
2j + 1 sites, 



^{x,t + l)= £ R^ mf ^(x + 2m,t), t = 0,l,2 



(j) „/,0')/ 



(1) 



We use Wigner's rotation matrices [16|, Il7| ] specified by three real parameters called the 



Euler angles a,/3, and 7 as the quantum-coin matrix, 



-j < m,m' < j 



with 



2j+m~m' -21 



sin 
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2e+m'-m 



Here 

yj (j + mjUj - mjUj + m')\(j - to')! 



T(j,m,m',£) 



' (j - ml - £)\(j + m- £)W.(£ + m'- to)!' 
and the summation Y^,£ extends over all integers £, for which the arguments of the factorials 
are positive or null (0! = 1). The position of the quantum walker at time t is denoted by 
Xf and the probability to find a walker at site x at time t is given by 

P {j \x,t) = Prob(X 4 (i) = x) = \^ j \x,t)\ 2 . 

The ratio X^/t is called the pseudo velocity of walker 13] and its r-th moment is given as 




xez 



at each time t = 0, 1, 2, • • •. 



For simplicity, we will assume that at the initial time t = one quantum walker exists at 
the origin, 

¥i\x,0)= ( j> { 3) 5 xfi 

with 

0o ) = T (lj, Qj-h Q-j+i, (2) 

where q m G C = the set of all complex numbers, — j < to < j, with J2m=-j \Qm\ 2 — 1- m 
this paper the left-superscript T denotes the transpose of vector or matrix. We usually call 
a (2j + l)-component qudit, which the quantum walker possesses at t — 0, 
For \a\ < 1, let 



M*; a) = -7j \ A o ^lxKlai}, (3) 

where l^} is the indicator function of a condition tu; = 1 if the condition u> is satisfied, 
= otherwise. We call it Konno's density function, since it was first introduced by 
Konno to describe the limit distributions of the standard two-component quantum walks in 
his weak limit-theorem fiol 
proved, 

dvv r i/W)(v), r = 0,1,2, 



ll\ . In an earlier paper [l4j . the following limit theorem was 

or 




with 



= E (^;cosf ) ^ (£) + i {w+1)is odd) A«^(,), (4) 



where A4^>' m '(x) are polynomials of x of order 2j. That is, the long-time limit distribution 
of pseudovelocity of quantum walker is described by superposition of appropriately scaled 
Konno's density functions ([3]) with weight functions /A^' m \v/2m), and a point mass at the 
origin with intensity if the number of states (2j + 1) is odd. In the previous paper jjjj], 
however, explicit expressions for the weight functions Ai^' m \x) are given only for j = 1/2, 1, 
and 3/2, since the functions seem to become very complicated as the value of j increases. 

In the present paper we will calculate Ai^' m \x) for large values of j and study the 
asymptotics of the limit distributions (HI) in the j — > oo limit. We will report our observa- 
tion that, when j becomes very large, complicated structures of the limit distributions are 
smeared out and a universal monotonic convex-structure appears around the origin. This 
observation leads us to a discussion on a possible route from quantum walks to classical dif- 
fusion. Relationship between the quantum-walk behavior and diffusive behavior of classical 



random- walk is an important topic in the study of quantum walks [181 ]. 

This paper is organized as follows. In Sec. II, hermitian matrix-representations of weight 
functions 'M.^' m \x) are introduced and formulas are given for the matrix elements, which 
are useful to calculate the limit density-functions for large values of j. We analyze limit 
density-functions for large j in Sec. Ill and clarify the properties of convex structure, which 
appears around the origin in limit distributions for sufficiently large values of j. Crossover 
phenomenon from quantum walks to classical diffusion associated with the j ^ oo limit is 
discussed in Sec. IV. Appendices are used for some details of calculations. 

II. WEIGHT FUNCTIONS 

A. Hermitian-matrix representations 

We note that the weight functions Ai^' m \x) are represented using (2j + 1) x (2j + 1) 
hermitian matrices M^' m ^(x) and (2j + l)-component initial-qudit (j2J) as 

M {hm \x) = [^[M^ix)}^ . 

In Appendix A, we give the matrices M^"%) = (M^ 2 (x)) for j = 1/2, 1,3/2 as exam- 
ples. In general, 

M mi =M ii -3<mi,m2<j, (hermitian condition) (5) 
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and 

M { l^_ mi (x) = (-l) mi ^ +2m MS:^(-x), -J < m x ,m 2 < j. (6) 
We found that, if the indices mi and m 2 satisfy the condition 

mi < m 2 and vn\ > —rri2, (7) 

we have the expression 

2 <i < 2 



.O'.m.mj) R (m,m 2 ) 

*l.<2 *i.*2 //l(i. m . m i)\ /n(m,m2)^ 



X 

fc 1= fc 2 =0 



Here the summations an d X^ 2 extend over all integers of t\ and £ 2 , for which the 
arguments of the factorials are positive or null, 

and 

= E E E Ll ) (1°) (I 1 ) (-l)*°+^ 2 ^)*«- 2 ^) (9) 

fco=0 fci=0 fc 2 =0 \ 0/ \ 1/ \ 2/ 

with r = tan(/5/2), where [z] denotes the integer not greater than z. For example, f^(x) = 
rx, fW( x ) = (2r 2 + l)x 2 - 1, ff\x) = (4r 3 + 3r)x 3 - 3rx, ff\x) = (8r 4 + 8r 2 + l)x 4 - 
(8r 2 + 2)x 2 + 1. Note that f^ a \x) is even (respectively, odd) if a is even (respectively, odd' . 



The derivation of ([8]) is tedious but straightforward following the method given in [13j, [14| . 
The key formulas are found in Appendix C of [14J]. Combination of the expression ([8]) with 
the symmetry properties ([5]) and determines all elements of the matrix M^ ,w ^ (x) for any 
given j e {1/2, 1,3/2, ■■■} and m e {-j, + ■■■,]}. 

B. Recurrence formulas 

As shown in Appendix B, from our expression ([8]), we can derive the following recurrence 
formula for matrix elements M^^x), when the condition (CO) is satisfied, 

+^f? ] \^ 2in M m ,=-^h (10) 



where 

2j 



c(j;m 1 ,m 2 ) 



\J(j + m 1 )(J + m 2 ) 

By solving this recurrence formula under the initial condition flAlfl . and by using the sym- 
metry properties ([5]) and ([6]), matrices My*" (a;) can be easily calculated even for large j. 
Moreover, ([8]) gives the following relation, 

M^-V(x) = 2(jX + mi)(jX + m2) M^) (x) (11) 
which enables us to determine M.^'^ 1 ' (x) from 'M.^'^(x). 



III. ANALYSIS OF LIMIT DISTRIBUTIONS FOR LARGE j 
A. Numerical evaluation of exact formula 

Though the formula (JE]) with ([9]) is rather complicated, it is exact for any given values 
of j, m, nil and m 2 . Therefore, if we fix the value of x, it is easy to evaluate M^^, (x) 
with any precision using computer. In order to demonstrate validity of this procedure, here 
we compare the numerical evaluations of our exact formula ([8]) with the results of direct 



computer-simulations 14[ of the quantum- walk models for a relatively large value of j. As 
an example, we set j = 11/2 and (a, 0, 7) = (0,7r/2,7r). The initial qudit has 2j + 1 = 12 
complex components. Figure 1(a) shows the result, when we choose the initial qudit as 
(f) = T ((l + i)/2, 0,0, 0,0, 0,0, 0,0, 0,0, (1 - i)/2), where the thick lines show the exact 
limit-distribution of pseudovelocity obtained by the above mentioned numerical 

calculation and the scattering dots indicate the distribution of X t /t at time step t = 100 
obtained by direct computer-simulation. For this initial qudit, the limit probability-density 
j/C 11 / 2 )^) i s symmetric and well describes the distribution of X t /t for large t. If we choose 
the initial qudit as (p = T (l + i, 0, 1 + i, 1, i, i, 1 + i, i, i, 1 + i, i, 1 + i) the distribution becomes 
asymmetric as shown by Fig. 1(b). 

Figure 1(a) shows that the limit distribution for the j = 11/2 model (twelve-component 
model) is given by superposition of six Konno's density functions ([3]), each of which is 
appropriately scaled according to Eq.(Tj0). In addition to them, inside of the innermost 
Konno's density function, we can see a convex structure in Fig. 1(a) 
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FIG. 1: (Color online) Comparison between direct simulation results and the probability densities 
of limit distributions for the twelve-component model, (a) Symmetric and (b) asymmetric cases. 

B. The convex structure 

From now on we fix the parameters of quantum coins as a = 7 = and the form of 
the initial qudits as 0o = T (q, 0, • • • , 0, q) with q = (1 + i)/2. In order to avoid Dirac's 
delta-function peaks, we will assume that the number of states 2j + 1 is even. Now we see 
the j-dependence of central convex structures. Figure 2 shows the central parts of limit 
distributions with a fixed window, X t /t G [—2,2], for variety of j's, when we set (3 = 7r/2. 
When 2j + 1 > 8, we see convex structures for (3 = tt/2. 



to 




FIG. 2: (Color online) Central parts of limit distributions for the models with 2j + 1 = 6 (plotted 
by crosses), 8 (triangles), 14 (squares), and 50 (circles). The convex structure appears around the 
origin, when the number of states with 2j + 1 becomes greater than eight for (3 = tt/2. 
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In order to verify the fact that, for any value of parameter /3, the convex-structure appears 
around the origin, if the number of state 2j + 1 becomes sufficiently large, we calculate the 
second derivative of the density of limit distribution (j3J). By using the exact expression of 
weight functions given in Subsection II. A, we have obtained the result, 



dv 2 



v=0 



0<m<j 



5m 



- cos 2 (/3/2) 

7TCOs(/?/2) 

1 



cos 2 (/3/2) 



2(2m 2 -j; 



(2j)! 



2 2 i- 1 (j + m)\(j -m)\ 



(12) 
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FIG. 3: (Color online) Values of second derivatives of limit density-functions at the origin are 
plotted for various f3. Crosses are for = vr/10, squares for (5 = ir/2, circles for (3 = 227r/25, and 
triangles for (3 = 467r/50, respectively. 



In Fig. 3 we plot the values of (fT2"j) with changing values of j for (3 = vr/10, tt/2, 227r/25 
and 477r/50. As demonstrated by this figure, we can prove that for any given (3 G [0, 7r], 
these is a critical value j c (P) such that d 2 v^\v)/dv 2 \ v =o < for all j > j c (P)- 



C. Smoothing by weight function in large j 

As given by Eq.(j4]), the probability density of limit distribution of pseudovelocities is given 
by superposing Konno's density functions ([3]) appropriately scaled. Since each Konno's den- 
sity function /i(t>/2m; cos(/3/2)) has a finite open support v G (— 2m cos(/9/2), 2m cos(/?/2)) 
and diverges at the boundaries of it, we see pikes at v = ±2m cos(/?/2), < m < j, in 



S 



the limit distribution as shown by Fig. 1(a), which was given for the case 2j + 1 = 12 and 
P = tt/2. 

For a given value of parameter (5, however, if we set j 3> j c (P), limit distributions seem 
to be quite different from that shown in Fig. 1(a). Figure 4(a) shows the limit distribution 
for the 2j + l — 50 case with (3 = tt/2. Pikes vanish in a central region and the convex 
structure at the origin becomes very evident. 
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FIG. 4: (Color online) (a) Probability density of limit distribution for the case 2j + l = 50, (5 = tt/2. 
(b) Functions H w / 2 ( m ) plotted for the cases 2j + 1 = 50 (crosses), 96 (squares), and 130 (circles). 



This smoothing phenomenon is caused by the interesting property of weight func- 
tions AiV' m >; for large j, A4^' m \v/2m) becomes to attain zeros at the points v = 
±2m cos(/5/2), < m < j, where the scaled Konno's density functions, \x{vj2m\ cos(/?/2)), 
diverge. In order to see this fact, for given j and j3, we define a function of m by 

Hf(m) = M {l ' m) {cos(3/2) 

(2j)!^^ (-l)^(cos(/3/2)) fc ^ 
" 2^~ x ^Q^o h\(j + m-k 1 )\k 2 \(j -m-k 2 )\ even} > 

< m < j. Figure 4(b) shows the functions, when (3 = tt/2, for the cases 2j + 1 = 50,96 
and 130. As increasing j, the central region, where Hg\m) = 0, becomes wider. 
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D. Rescaling of limit density-function 



By definition of our models, the range of elementary hopping of quantum walker at each 
time step is 2j; see Eq.(jT|. Then distribution of pseudovelocities of quantum walker spreads 
in an interval (— 2j cos(/3/2), 2j cos(/?/2)). 

In order to discuss the j — > oo limit of the series of our models, here we introduce the 
rescaled variable 



X, 



U) 



0') 



2jcos(/?/2) (13) 
for each value of j3. Figure 5(a) shows the limit density-functions of the rescaled pseudoveloc- 
ities Xf /t for 2j + 1 = 10, 20 and 50. In this variable, the support of the limit distribution 
is fixed to be (—1, 1). As shown by Fig. 5(a), the central convex structure becomes sharper 
monotonically in increasing the value of j. Corresponding to (fl~3l) . we plot (jjHp^s as func- 
tions of m/(jj, where Uj = y/2j, in Fig. 5(b) for (5 — ir/2. It is interesting to see that the 
locations, where the functions take non-zero values, are now fixed. 
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FIG. 5: (Color online) (a) Central parts of density functions of limit distributions of rescaled 
pseudovelocities for 2 j + 1 = 10 (crosses), 20 (squares), and 50 (circles), (b) OjHp(m) versus m/aj 
with ex,- = y/2j for 2 j + 1 = 50 (crosses), 96 (squares) and 130 (circles), when = tt/2. 



IV. DISCUSSION 

Now we discuss the j ; — > oo limit of the present series of one-dimensional quantum-walk 
models. If we consider the situation that j is finite but very large, the limit distribution 
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of the rescaled pseudovelocity X^ /t will have a simple profile. There will be a sharp 
convex structure at the origin and accumulations of pikes only in the very vicinity of the 
boundary points -1 and 1 of the support of limit density-function, but in other regions in 
the support the density of distribution will be almost zero. In the j — > oo limit, the central 
convex-structure will become a single point mass (i.e. a Dirac delta-function) at the origin. 
It implies that the pseudovelocity is almost zero, that is, the quantum walk will lose its 
velocity in the j — > oo limit. 

In the series of one-dimensional quantum-walk models studied in the present paper, we 
have used Wigner's rotation matrices Ry' as quantum coins. We should note that is 



introduced 



16 



171 ]. when rotations in the three-dimensional real space are quantized and 



the index j is defined as a quantum number, which specify physical states allowed in the 
quantum mechanics. 

In the standard quantum mechanics, a small but finite parameter h (the Planck constant 
divided by 2n) is introduced and physical quantities are quantized to have only discrete 
values of the form fi x (quantum numbers). For example the energy levels of a harmonic 
oscillator with the angular frequency uj are given by E n = hu(n + 1/2), n = 0, 1, 2, . . ., and 
the square of angular momentum L and its ^-component are given by L 2 = h 2 £(£ + 1) and 
L z = km, £ = 0, 1/2, 1, 3/2, • • ■ , m = —£, — £ + 1, ■ ■ • , £. If the quantum system involves Ti 
explicitly, it is easy to consider its classical correspondence by taking the limit h — > 0. Even 
though the system does not include the parameter h explicitly, the classical limit can be 
realized by taking a large quantum- number limit, since physical quantities should be given 
of the form % x (quantum numbers) as the above examples show. 

Then we can expect that, if we take j oo limit appropriately, the classical diffusive 
behavior will be observed [181 ] . One possibility is to see a crossover phenomenon from 
quantum-walk behavior to classical diffusion in j — > oo. Assume the form X^ ~ jtF(t/j e ) 
with a scaling function F(z) such that F(z) ~ z^ 1 ^ 2 in z — > and F(z) — > const, in z — > oo, 
where 9 is an exponent. For finite j, X^ ~ jt in t — * oo, as we have shown in the present 
paper. On the other hand, for finite t, we will have X t ~ x / 2 = j 1+6 ^ 2 y/i in 

j — > oo; that is, x\^ j ' j l+0 l 2 is diffusive in large t. It will be an interesting future problem to 
clarify the phenomena, which are realized when we take a proper classical limit in the series 
of quantum-walk models. 
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APPENDIX A: THE MATRICES M^ m )(a;) FOR j = 1/2,1,3/2 



Let t = tan(/3/2) and /«(x) = rx, fW(x) = (2r 2 + l)x 2 - 1. Then 



M (1/2 ' 1/2) i 



x) 



x) 



1 — x rxe 11 
Txe~ %1 1 + x 



1 1(1 -*) 2 

V2 



V2 



(l-x)fW(x)^ \f^\x)e 2i ^ 



(Al) 



^(l-x)/«(x)e- i7 (l-x)(l + x) 



\\tf\x)e-^ 



V2 
2 



(l + x)fW(x)e-*y Hl + xf 



J 



^-{l-x) 2 ft\x)e^ ^-{l-x)f?\x)e 2 ^ jf^\x)e 3 ^ 



\ 



^(l~ x ) 2 f^(x)e-^ |(1 -x) 2 (l + x) l(l-x)(l+x)fi 1 \x)e i ^ ^(1 + x)f?\x)e 2 ^ 

^{l-x)f?\x)e- 2 ^ l(l~x){l+x)f^\x)e^ |(l-x)(l+x) 2 ^(l+^/^We^ 



(l+x)fi 2) (x)e~ 2 ^ 



^(l + x) 2 fi 1 \x)e-^ 1(1 + x) 3 



/ 



APPENDIX B: DERIVATION OF EQ.fllO]) 



When m x ^ —j, Eq.(jS]) gives 



1 1 \ / 1 11 

xV[j--,m 1 --,j--,0)v(j-- 







x 



''g^/j+mi lUj 7 ^V_ 1 )fc 1(a .fci+* a _ a .fc 1 +* a +i). 
fel =o fc 2 =o V fc i J\k 2 J 

The summations over &i and &2 are carried out 



as 



j—m,2 

E 

fc 2 =o 



fci=0 



E 

fci=0 



(-l) fcl x fcl+1 
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E ,i Ep;r i )( , '-;j(-i) t '^ +fe 



fe 1= fc 2 =0 

We note the equality 



, / 1 1 1 \ / 1 1 1 \ 

c{3;m u m 2 )T l^j -~,m x - -,j - -, OJ T \ j - -, m 2 - -,j - ~,0J 

= r(j,mi, j,o)r(j,m 2 , i,o). 



Then Eg. (1101) is derived. When m x = — j, it is enough to only consider the case mi = j 
under the condition ©. This case is trivial; M^\x) = 2~ 2j+1 f^ 2j \x)e' 2ij ' / . 
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